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Abstract 

We apply the reduced phase space quantization to the Kasner universe. We construct the kine- 
matical phase space, find solutions to the Hamilton equations of motion, identify Dirac observables 
and arrive at physical solutions in terms of Dirac observables and an internal clock. We obtain 
the physical Hilbert space, which is the carrier space of the self-adjoint representation of the Dirac 
observables. Then we discuss the problem of time. We demonstrate that the inclusion of evolu- 
tion in a gravitational system, at classical level as well as at quantum level, leads respectively to 
canonically and unitarily inequivalent theories. The example of Hubble operator in two different 
clock variables and with two distinct spectra is given. 
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I. INTRODUCTION 



The standard model of cosmology is based on the remarkably simple solution to general 
relativity, the Friedman-Robertson- Walker (FRW) universe. It is expected, however, that 
a slightly perturbed FRW universe, when evolved back in time, at some moment close 
enough to the big bang singularity, will lose its space-like symmetries. Therefore, in order to 
understand the singular conditions from which the universe emerged nearly 14 billion years 
ago, a study of more general cosmological spacetimes is needed. 

A general solution of general relativity in the vicinity of cosmological singularity has 
been studied by Belinskii, Khalatnikov and Lifshitz (BKL) in In the BKL scenario, as 
spacetime approaches singularity, the time derivatives of the gravitational field are shown to 
dominate over all spatial derivatives for relatively long stretches of time. Surprisingly, the 
evolution of the general gravitational field turns out to be well approximated, at each point 
separately, by a sequence of the so-called Kasner epochs. Each epoch is a vacuum solution to 
the homogenous spacetime model of Bianchi I type. The transitions between epochs are the 
effect of non-negligible spatial curvature, which arises quickly and vanishes after a relatively 
short period of time. In the BKL scenario the universe undergoes an infinite number of 
chaotic-like transitions and eventually collapses into a singularity in a finite proper time. 

It is commonly believed that the incompleteness of classical theory, which breaks down 
at the singularity, will be overcome by quantization of the gravitational degrees of freedom. 
For this purpose, the Dirac method of quantization is usually employed (see e.g. [2-01 )• 
In this paper we focus on quantum theory of the Kasner universe. We follow, however, 
an alternative way to quantum theory, namely the reduced phase space quantization (see 
e.g. [HI, 0]). Since the Kasner model plays a central role in the BKL description of a 
generic cosmological singularity, we believe that the present and future investigations into 
this model supported by current and forthcoming astrophysical and cosmological data can 
help obtaining new insights into the universe's origin. 

The fact that in canonical general relativity the evolution of gravitational fields coincides 
with gauge transformation or equivalently that there is no privileged time standard to mea- 
sure motion, gives rise to the so-called problem of time (see e.g. Kuchar [3]). The problem 
of time consists of a few related though distinguishable issues. Following Kuchar's termi- 
nology, we will treat in this paper the most fundamental issue, namely the multiple choice 
problem. In essence, it states that two different choices of time may produce different quan- 
tum theories. We will show how severe the problem is and examine its origin. In the view of 
the results obtained in this paper, the proposals for explaining the multiple choice problem 
existing in the literature are unsatisfactory. For example, in Isham ^ we can read that 
two different choices of time lead to canonically equivalent theories, which admit unitarily 
inequivalent quantum representations due to the Van-Hove phenomenon joj. We will show 
that in fact the multiple choice problem can be traced back to the canonical formulation of 
general relativity and thus studied at classical level. 

The application of the reduced phase space approach to the Kasner universe turns out 
to be manageable and quite straightforward. The system consists of a single constraint on 
the kinematical phase space, which is six- dimensional. We identify the (physical) reduced 
phase space, that is the space of Dirac observables, which is four-dimensional. We consider 
two examples of clock variables to introduce the physical evolution and construct the so- 
called true Hamiltonians. We argue that the evolution is free of the singularity present in 
the classical theory. We compare the spectra of the Hubble operator in two different clock 
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variables and show that they are very different. 

Before we start let us introduce the notation that will be used throughout the text. 
The canonical variables, which follow from the Legendre mapping applied to the Einstein- 
Hilbert action, parametrize the kinematical phase space, denoted by V. In this space, 
the Hamiltonian constraint H is introduced. The constraint surface, defined by = 0, is 
denoted by S. The Dirac observables, denoted by Pj, are defined by the relation {Dj, H} = 
and with the domain restricted to the constraint surface. The space of Dirac observables, 
called the reduced phase space, or the physical phase space, will be denoted by Vr. In 
Rovelli it is proposed to call kinematical phase space functions by partial observables 
in order to emphasize that they can be measured by observers, though the outcome of 
such a measurement made alone cannot be predicted by theory. There is, however, one 
kinematical degree of freedom for which the theory makes the immediate prediction, i.e. the 
Hamiltonian constraint vanishes, H = 0. In what follows, by partial observables we mean 
functions restricted to S (which seems to be slightly different from Rovelli's notion), and 
which will be denoted by Vi. 



II. LAGRANGIAN FORMULATION 



The Hilbert-Einstein action reads 



She = I RV~9 d'^xdt 
'ncM 



(1) 



where g, R are the metric determinant and the Ricci scalar, respectively. The integral is 
taken over an open subset Q of the manifold Ai. We specify action ([1]) to the case of 
vacuum Bianchi I model with = i? x E, where E is a compact spacelike leaf and we 
assume the following metric type: 



ds'^ = -N\t)dt^ + al{t){dx^f + al{t){dx^f + al{t){dxy 
which leads ([1]) to the form (see Appendix 



She = 2 / dt I Naia2a^ 
'to Jt. 



^ ATnA ATI ^ ^ 



ai 
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Applying the variational principle to (|3]) gives the Lagrange equations: 



d^x 



(2) 



(3) 



ft ^ 



(4) 



The solutions were first found by Kasner in Later, they were rediscovered by Taub in 
12| . who gave the solutions in the following form: 



ds^ = -d^ + e^'dx^ + t^P^dy"^ + t^PHz^ 
where the constants pi, p2 and ps satisfy: 



(5) 



(6) 
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All the above solutions, except for = 1, p2 = 1 or = 1, admit a cosmological singularity 
for t = 0. 

It is easily seen that action ([3]), equations of motion (j4]) and solution ([5]) are invariant 
under any time re-parameterization such that t i{t) and i— )■ ^N. This is the only 
gauge freedom, which is preserved under the reduction of the Hilbert-Einstein action to the 
homogenous spacetime. Therefore, by gauge transformation we mean any change of time 
parameter t ^ t{t). 



III. HAMILTONIAN FORMULATION 

Using action Q we define the momenta0: 

N 9L ^ dL 2 ■ 

The Legendre mapping ([7j) is singular, and its range is the submanifold of the phase space 
given by vr^ = 0. The Hamiltonian reads: 



The Dirac analysis [3] leads to the reduction of phase space by the conjugate pair (A^, 
and the introduction of the Hamiltonian constraint: 



AT i 

if := - V ^ (ttV, - TT^a, - Tr'^a,) ^ (9) 

where now the lapse is a Lagrange multiplier. One may verify that the vanishing of 
Hamiltonian constraint (|9]) and Hamilton's equations are equivalent to the first and the 
second of the Euler-Lagrange equations in (jl]), respectively. 



A. New variables and motion 



We introduce the new canonical variables: 

X, := ^a.afc , P, := -4^ (10) 
in which the Hamiltonian constraint dH]) reads: 

The symplectic form on the kinematical phase space is defined as: 

u■=^dX,^ dPi (12) 



For the sake of simplicity, from now on we drop the integration over the compact space-like leaf S whenever 
it should appear. 



4 



and its minus inverse is the Poisson bracket, i.e.: 



i = \-lJ_l U_) (13) 



The Hamilton equations in the gauge = A^/2^/X[X2X^ read: 

^^=QP^=-UX,P, + X,P,) (14) 

Pi = = n{x,Pj + XkPk) (15) 

From combining the above equations into one: 

XiP, + XiPi = {X,P,),t = (16) 

we obtain that Fj := XjPj are constants of motion (as we will see in a moment they may be 
identified with some of the Dirac observables) . Putting this back to (fT4l) . f|T5l) we find that: 

X,{t) = X,(to)e-(r^+r^)(*-*°) , P.(t) = P,(to)e(r^+r'=)(*-*«) (17) 

where Pi{to) and Xi{to) are the initial conditions for the Hamilton equations and to "will be 
specified later. The physical solutions (|T71) should satisfy the Hamiltonian constraint (fTTil . 
which can be rewritten now as: 

5^F,F,^0 (18) 

i>j 

The change of the arrow of time t i— )■ — t in f|T7|) is equivalent to the sign change Fj i— ?■ — Fj 
for all z's, so we can add the condition: 



5^ F. > (19) 



which ensures that the singularity is approached as the time t grows. Moreover, the following 
three cases: 

Fi = F2 = 0, F2 = F3 = 0, Fi = F3 = (20) 

can be shown to correspond to the Milne space, which can be isometrically embedded in 
Minkowski spacetime and thus are non-singular (i.e. the coordinates are singular, not the 
spacetime itself). We exclude them from the phase space. 



B. Dirac's observables 

It is known that for the Hamiltonian satisfying dH 7^ in a neighborhood of the constraint 
surface H = 0, one may locally introduce such a canonical parametrization of the kinematical 
phase space that the canonical coordinates: 

(X„F,), 1 = 1,. ..,n (21) 

are replaced with the new canonical pairs: 

{H,T),{Xi,P,), i = l,...,n-l (22) 
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such that the variable T is canonically conjugate to H and the symplectic form now reads: 

u = dTAdH + dX^A dPi (23) 

It is now easily seen that the space of functions which commute with the Hamiltonian H is 
given by: 

{H,V,} = 0^V, = V,iX„P,,H) (24) 

which restricted to the constraint surface H = can be identified with the following space 
of functions: 

Vi = Vi{X„P,) (25) 

An easy way to find all the Dirac observables is by pulling back the symplectic form u to 
the constraint surface: 

u\h=o = dXi A dPi (26) 

and ensuring the pulled-back two-form is in canonical form. Obviously, this recipe does not 
depend on gauge as 

uj = (NH) -dTAd^ + ^-dTA d{NH) + dXi A dPi (27) 

so uj\nh=o = ijj\h=o for any 7^ 0. 

In what follows, we will obtain the complete set of Dirac observables by restricting the 
symplectic form oj introduced in (fT2|) to the constraint surface (fTSjl . Applying the mapping: 

1 1 

Ti = -X V + z 

2 2^ 

r2 = ^2: - - ^ (28) 
^35 

one arrives at the following form of the constraint f|T8l) : 



The constraint is solved 


2 2 2 n 
X — y — z =0 

in the new coordinates (r, </>) G i? x 5 such that: 


(29) 


so that Fj's read: 


X = r , y = r cos , z = r sin (p 


(30) 




Ti = 2'"(-'- ~ '^os + 2 sin 0) 


(31) 




r2 = 2^{^ ~ cos — 2 sin 0) 


(32) 


One checks that setting 


T3 = -^r^S + 5 cos 0) 
the arrow of time in fflQl) gives 


(33) 



5^r, = ^(7 + cos0)>O (34) 
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and leads to r > 0. The exclusion of the Milne space cases ( 120|) . which are equivalent to 
(/) = and cos = — |, restricts the parameter (j) in the following way: 

(O,0i)U((/)i,02)U(02,2^) (35) 

where (f)i and (p2 are the solutions to cos^ = — |. The starting kinematical phase space was 
parameterized by the six coordinates (Xj,Pj). Since Xi > 0, we could use the coordinates 
(Xj, Fj) as well. Thus, constraint surface flTSl) may be parameterized by the five coordinates 

The complete set of Dirac observables and their commutation relations can be found 
by restricting the symplectic form a; in "P to the constraint surface S. Let us denote the 
embedding of the constraint surface in the kinematical phase space hj E : S ^ V so that 
E* : C°°{V) 3 {P\ Xi) H- (r, 0, Xi) e C°°{S) is the restriction of the kinematical phase space 
functions to the constraint surface functions: 



E*{J2dX.AdP.)=E*{J2dX.A^) 

= d[{- — - cos + sin 0) InXi + (- — - cos — sin 0) lnX2 + (- + - cos (j)) InXs] A dr 

115 

+ (i[(-rsin0 + rcos(/)) InXi + (-rsin0 — rcos0) lnX2 rsin^lnXs] A (i0 

2 2 ^ 

= dOx Adr + d02 A d(j) (36) 

where we have defined: 

11 11 3 5 

Oi= (- - - COS + sin (j)) In Xi + {-^ - cos - sin 0) In X2 + (t + T ^) ^3(37) 

115 
02= {-r sin + r cos 0) In Xi + (-r sin — r cos 0) In X2 — —r sin In X3 (38) 
2 2 ^ 

The space of Dirac's observables is called the reduced phase space Vr : (r, 0, Ci, C2) G 
X IiU I2U I3 X R X R. Apparently, the space of Dirac observables in the Kasner universe 
is not simply connected. The four Dirac observables together with any constraint surface 
function t such that 

{t,H}\^^^y^O (39) 

form a coordinate system {r,(j),0i,02,t) on the five-dimensional constraint surface S 
equipped with the two-form: 



Us 



:= dX, A dPi 



= dOi Adr + d02 A dcj). (40) 

H=0 



induced from the symplectic form u on the kinematical phase space V. It should be added 
that the freedom in the choice of the fifth coordinate t on iS is bigger then the freedom 
in choosing the lapse function N : S ^ which fixes only the first derivative of t with 
respect to the Hamiltonian vector field, i.e. {t, H}\^_^ = N~^, where H itself is taken with 
the lapse equal 1. 
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C. Partial observables 



Solutions found in (1171) include both the physical and non-physical sector. In the physical 
sector, however, all the solutions should be expressible in terms of coordinates on S. Setting 
-^i(i^o) = 1, we express all the constants occurring in ( |T7I) . that is Fj, Xi(to), Pi(to), in terms 
of the Dirac observables Oi,02,4>,r and arrive at: 



Pi 


= -T^i^ — COS + 2 sin ( 


^) exp(-r(5 + 3 cos — 4 sin 0) {t — 


to)) 


(41) 








50i sin 4>+ 1 O2 (3+cos q 




P2 


= — cos (p — 2 sin i 


i)) exp(-r(5 + 3 cos + 4 sin (j)){t — 


to^^G —5 — 3 cos <^+4sin cj) 


-(42) 




1 


O1 (2 sin 0-4 cos 0)4 


- (-2 + 2 cos tp+A sin <p) 


(43) 


P3 


= -'"(3 + 5 cos <p) exp(r(l — cos (p){t — tojje -s-^' 


cos (/)+4 sin (j) 




= exp(— -r(5 + 3 cos (p 


-4sin0)(t-to)) 




(44) 






50i sin 0+ i O2 (3+cos 


) 


(45) 


X2 


= exp(— -r(5 + 3 cos (/) 


+ 4sin0)(t - to))e -5-3 cos 0+4 sin 






Oi (2 sin 0-4 cos 0)+ i O2 (-2 + 2 cos 0+4 sin 


10) 


(46) 




= exp(— r(l — cos0)(t 


— i^Q^^C —5 — 3 cos 0+4 sin 





The above solutions are solutions (|T71) restricted to the constraint surface S. The variable 
t parameterizes the gauge orbits in the constraint surface S. After the specification of the 
value of to for each gauge orbit, that is defining to as a function of Dirac observables, the 
variable t becomes the fifth coordinate, which assigns a specific value to each point in S. 

In analogy to the Friedman cosmology, we will be interested in the Hubble and decelera- 
tion parameters, which in the Kasner model are introduced for each of the three directions 
(see Appendix and read: 

rfl — cos (/) + 2 sin (/)) fr \ °i 0-4 cos 0)+ 1 02(1+3 ■^+4 sin 0) 

Hi = — ^ ^ — exp ( - (7 + cos (j)){t-to)\ e io+6cos0-8sin0 (^47) 

r(l — cos (/) — 2 sin (/)) fr \ C'l(7sin0-4cos0)+iO2(l + 3t:oS'*+4sin'P) 

H2 = ^ ^ ^exp ( ^(7 + cos0)(t - to) 1 e io+6cos 0-8sin0 (^43) 

3 -|- 5 COS I T \ (7sin (ji — 4 COS 0) sin 0+iC>2 (1 + 3 cos sin 0) 

= — — -^exp ( -(7 + cos0)(t - to) je io+6cos0-8sin0 ^49) 

f- 7 + cos (h ^ 

^1 = ^ 1 jrf— A - 1 ^0 

1 — cos + 2 sm (f) 
/7T 7 + cos 

^2 = V 2 — -— — - 1 51) 

1 — cos (p — 2 sm 

r- 7+COS0 , , 

ga = 2V2 \ - 1 52 

3 + 5 cos 



We note that the deceleration parameters Qi = 2a/2(^;^) — 1 are constants of motion, 
i.e. Dirac's observables. All the above six quantities determine and are determined by the 
components of the connection and curvature matrices (see Appendix |A]). Therefore, they 
represent the local properties of the Kasner universe and do not form the complete space of 
observables in the compact universe. 
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D. Gauge transformations 



For t to be a function on S, we need to specify to as a function of Dirac observables. We 
note that the choice: 

Ci(7sin0 - 4cos0) + iC2(l + 3cos0 + 4sin0) 
^° |(7 + cos0)(5 + 3cos(/)-4sin0) ^^^^ 

simphfies nicely the formulae for the Hubble parameters: 

r(l — C0S(?!» + 2sino!)) fr ,, \ . 

Hi = ^ ^exp ( -(7 + cos0)tj (54) 

r(l — cos (^) — 2 sin 0) fr, , \ , , 

H2 = ^ ^exp l-{7 + cos<P)t] (55) 

r(3 + 5cos0) \ 
Hs = ^^-^ exp (-{7 + cos (f>)t] (56) 

and leaves the formulas for deceleration parameters unchanged. 

We want, however, more than that and apart from simplicity we require that the time 
coordinate satisfies two extra conditions: (a) has a clear physical meaning and (b) the 
singularity occurs at its finite value. A distinguished choice is the cosmological time, for 
which the lapse function Ncos — 1 and the singularity is reached at tcos — (for all gauge 
orbits). 

We use the relation Ndt = Ncosdtcos to obtain the formula: 

'dtcos /■* 



t^^ = / -Jft dt = / -— dt (57) 

where t = tg defines the four-dimensional boundary of the constraint surface S, at which 
the singularity occurs. Then we insert N — A:\f2^JX^X2Xl and obtain 

-16v^ ( r ,^ \ 

r[" + cos 0) \ 4 J 

This time redefinition simplifies the formulae for the Hubble parameters even further: 

= + (59) 

V2(7 + COS0) tcos 

V2(7 + COS0) tcos 
H, = ' + '^"^'^ (61) 

2v^(7 + COS<^) ieo. 

It should be stressed that tcos occurring in the above formulae is a function on S. Thus, the 
cosmological time tcos may be related to the kinematical phase space functions. It cannot be 
done uniquely as there is no natural projection from the kinematical space to the constraint 
surface. However, there exists a class of kinematical phase space functions, which coincide 
with tcos on the constrained surface. To sum up, tcos is a partial observable as much as any 
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IV. PROBLEM OF TIME 



All internal clocks are given by the formula: 

rtcos 

T= / N-^dt,os + ro (62) 
Jo 

where N{tcos, T^i) > is a function of the cosmological time and Dirac observables and To(I^j) 
is any function of Dirac observables. Due to the change of time parameter 1 1— >■ r, any time- 
dependent function on the constraint surface, i.e. any partial observable, is transformed 
accordingly: 

V{V,,t) ^ ViVi, t) = V{Vi, t{T, V,)) (63) 

Thus, in general, V and V will have different dependance on Dirac observables and after 
quantization they may have different spectra. Therefore, it is meaningless to speak about 
spectra of partial observables, like energy density, curvature or volume, without reference 
to the choice of internal clock. To support this statement, it is enough to note that the 
commutation relations: 

{Vl,V2} ^ {Vl,V2} = {Vl,V2} + -^{t{T,Vi),V2} + -^{V,,tiT,Vi)} (64) 

becomes altered after the change of time. Note that the commutation relations between the 
Dirac observables can be obtained from the form a;^ given in (HOl) . where any variable chosen 
to play a role of time is interpreted as an external parameter. Thus, in the constraint surface 
the gauge transformation is not canonical and consequently it cannot be unitary in quantum 
theory. This implicates that the spectral properties of the same operator in different gauges 
should be, in general, different. 

The gauge transformation has another interesting feature. Suppose that we are given any 
two partial observables, which monotonically increase with time and whose range is identical. 
For example, let it be the curvature C and energy density p. Then there is always such a 
gauge transformation t ^ r' that energy density in one gauge is functionally identical with 
the curvature in another gauge, that is p\r = C\r'- In effect, the spectral properties can be 
shared by many different partial observables. If we drop the assumption about the identity 
of ranges, the two different partial observables will still share the identical dependence on 
Dirac observables for some time during the evolution. 

Obviously, one is unable to study the cosmological singularity problem without the notion 
of evolution, for example just by the inspection of Dirac observables. Thus, this gives rise 
to the very interesting question of the possible dependance between the choice of time and 
the fate of singularity in quantum theory. 



A. Geometrical formulation of the problem of time 

Let us take a closer look at the structure of systems with a Hamiltonian constraint. The 
constraint equation H = defines the embedding E : S ^ V of the constraint surface S into 
the kinematical phase space V having its dimension increased by 1. The constraint surface 
S is therefore odd-dimensional. All the physical motion takes place in the surface and the 
vectors tangent to the trajectories (gauge orbits) are given by the Hamiltonian vector field 
Xh (gauge generator), defined in the standard way: 

Xh: u{-,XH) = -dH (65) 
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where u is the symplectic form in V and the Hamiltonian vector field Xh\^_q, which is 
restricted to the constraint surface, will be denoted by Xh for brevity. 

The constraint surface is equipped with the two-form us induced from the kinematical 
phase space, us '■= E*u. The form a;^ is a singular closed two- form of maximal rank. It 
is represented at each point p G 5 by an antisymmetric matrix us{p) '■ TpS i— )■ T*S in the 
tangent space of the odd- dimensional manifold. Its null vector is Xh, satisfying u}s{-, Xh) = 
0. The null vector Xh is a generator of the line bundle with the projection (submersion): 
TT : I— )■ Vr from the constraint surface to the reduced phase space. The puUback tt* is 
understood as an injection from the space of functions Vi such that XH{T>i) = to the space 
of all observables. The diagram below illustrates the relations between V, S and Vr- 



_ Embedding, E ^ Projection, tt _ /r^r'S 

V i S > Vr (66) 

Since the form us is singular it cannot be inverted into the Poisson structure (which 
should be quantized). Suppose we introduce a projection on the tangent space at p E S, 
denoted by Pp : TpS TpS, such that Pp{X) = <S=^ X ~ Xh{p)- Then, the matrix ujs{p) 
can be inverted on the restricted domain {X G TpS : Pp{X) = X}. Though it enables to 
derive the inverse of us, this construction is ambiguous due to ambiguity in the choice of 
the projection Pp. Now, the (minus) inverse —ojg^p '■ TpS h- >• TpS is the Poisson bracket. 

The generalization of the above construction follows straightforwardly. Suppose we attach 
to each point g G iS a projection Pg of the considered type. In addition, we assume that 
there exists a slicing of S such that the tangent space to a slice at any q E S is identical 
with the range of Pg. Let the slicing be given by a function t : S ^ R. Now the Poisson 
structure for each q E S is given by the inverse of the form Us restricted to the hypersurface 
t = const containing q, and let it be denoted by —oo^]. 

The construction of —ujgl may be achieved by exploiting the simple relation that deter- 
mines the induced Poisson structure: 



-Uslit,Vi):={t,Vi}s,t = (67) 

with all other commutation relations being fixed uniquely by the independent-of-the-choice- 
of-time commutation relation between Dirac observables. 



B. Canonical transformations in constraint surface 



Once a slicing is introduced in the constraint surface, we arrive at the triple {S,uJs,t). 
Such a structure is well known in classical mechanics and is called a contact manifold. 
The restriction of us to the constant time hypersurfaces forms a symplectic submanifold 
with the symplectic form us,t- The form us,t can be pulled back to the constraint surface, 
LOs,t = t*us,t- It should be noted that in general us,t 7^ i^s- The relation between us,t and 
Us is the subject of the theory of canonical transformations and will be discussed below. 

One may think of a canonical transformation as a symplectomorphism, however, the 
contact manifold provides a better avenue to define this notion. Let us cite the definition of 
the canonical transformation from Abraham and Marsden [l^ . 

Definition. Let (Vi,ui) and (^2,1^2) be symplectic manifolds and {R x Vi,Cji) the 
corresponding contact manifolds. A smooth mapping F : R x Vi ^ R x V2 is called a 
canonical transformation if each of the following holds: 
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Dl: F is a diffeomorphism; 



D2: F preserves time; that is, F*t = t; 

D3: there is a function Kp £ C°°{R x Vi) such that F*U2 = u^k, where uk = Cji + dKp A dt. 

In what follows, we will use the symbol Ht instead of Kp- There are a few observations 
that can be made in effort to understand the above definition. First, we note that the slicing 
of the constraint surface: 

S (68) 

is needed in order to introduce the canonical transformations. However, the lack of this 
slicing is the essence of gauge invariance and it can only be postulated. Once it is done, any 
canonical transformation, according to condition D2, preserves chosen time. 

This leads us to the next observation that different choices of time must produce canon- 
ically inequivalent theories. It confirms our earlier result that the choice of slicing fixes the 
Poisson structure, which is related to the defining of an ambiguous procedure by which a 
non-invertible matrix (w^) can be inverted (w^]:). 

Another observation is as follows. Suppose there is a given slicing, t. The constraint 
surface can be parametrized as S = t x Vr, where Vr is the reduced phase space. We note 
that us^t = and that us = ujr{= t*UR). In this contact manifold there is no Hamiltonian, 
since the reduced phase space consists of Dirac observables for which = 0. Now we may 
introduce the evolution into the system by considering time-dependent reparametrization 
of the constant time hypersurfaces. As we will see below this will render a non-vanishing 
Hamiltonian. 

Consider the following canonical transformation F : S ^ S, such that F* : "Dj i— )■ 
Vi{q,p,t), F*t = t and: 

F*us = ujs,t + dHr A dt, (69) 

where u)s,t = i*^s,t is the puUback of the symplectic from ujs,t living in the leaf t = const, 
and parametrized with new coordinates {q,p). The new coordinates are in general time- 
dependent. 

In the coordinates [T>i,t), the null vector of the form a;^ is given by Xfj = dt. The 
canonical transformation F changes the coordinates and the coordinate expression for Xh 
accordingly: 



F*Xh = dt + XT = dt- us]{-, Ht) = dt + {-, Ht} 



(70) 



Now it is seen that the evolution of the observables in new coordinates is not only given 
through the explicit dependence on time t but also through the true Hamiltonian Ht. 

However, not all time-preserving diffeomorphisms are canonical transformations. There- 
fore it is useful to introduce an alternative formulation which relies on a generating function 
W: 

Pidq' + HTdt - F*eR = dW (71) 

where 6r, satisfying —dOR = ur, is the canonical (Poincare) one-form in the reduced phase 
space, and {pi, g*) are the new canonical pairs. is a generating function such thatH: 

dW dW 

W = W{q^^,q.,t), PD. = -g-r, P^ = ^ (72) 



We denote by q]^ and pn^i the reduced phase space basic variables. 
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and the relation between W and Ht reads: 

The final remark is that the true Hamiltonian is quite arbitrary and it does not depend 
on the particular choice of time but rather on the choice of basic variables on the constant 
time submanifolds once time is given. 

Let us sum up. In the kinematical phase space V, the symplectic form, u, and the Pois- 
son bracket, —uj~^, can be considered interchangeably. However, in the constraint surface 
S the induced two- form, wg, is singular and one cannot define the Poisson bracket in S un- 
ambiguously. The bracket is needed to compute the commutation relations between partial 
observables prior to quantization. The choice of the clock variable t : S ^ R, is equivalent 
to the choice of the following commutation relation: 

{t,Vi}s = (74) 

The last relation fixes the Poisson bracket between any pair of partial observables on S. 
In principle, one can choose any function t : S ^ R admitting Xnit) > to slice the 
constraint manifold. The fact, that different slicings lead to canonically inequivalent theories 
is confirmed by the formula fTMj) . 

It should be emphasized that the Poisson bracket between any pair of Dirac observables is 
given uniquely in a constrained system. There is the unique two-form in the reduced phase 
space, ur, such that its pullback to the constraint surface gives 'k*{ur) = us- The form ojji 
is invertible, since the reduced phase space Vr is even-dimensional, and the Poisson bracket 
can be computed. This is in agreement with the fact that the Poisson bracket between Dirac 
observables in the kinematical phase space, V, is given uniquely and independently of the 
choice of time. 



C. Problem of time in Dirac quantization 

In the Dirac quantization, the kinematical phase space is quantized so that the kinemat- 
ical Hilbert space is obtained and the Hamiltonian constraint is promoted to a self-adjoint 
operator H ^ H . Then the theory is constructed via the solutions to the quantum constraint 
equation, i.e. 

H^ = Q (75) 

The solutions il) normally do not belong to the kinematical Hilbert space, and the Hilbert 
space structure needs to be reintroduced. There is the idea, called 'deparameterization', 
which cures the problem and at the same time nicely introduces the concept of evolution of 
the system. The idea is to reformulate the constraint equation H = Q m such a way that the 



quantum constraint equation (175!) gets a Schrodinger-like form [ij, ll5| . This procedure is 



very closely related to another procedure used for the Dirac quantization, namely the group 



avaraging method [16[. Therefore, we will focus here only on the idea of deparametrization 
while keeping in mind that the problem of time is in fact method-independent and can be 
also formulated in the context of group averaging. 

Suppose that we consider a gravitational system including a scalar field. In this case (see 
e.g. [l5|), H K, p^^ — Cgr and we obtain a Schrodinger-like equation: 

- ihj-i^ = VdcR^ (76) 
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so that the scalar field, (f), plays a role of a time parameter. The non-vanishing true Hamil- 

tonian, ^/C^gri is expressed in terms of the rest of the kinematical degrees of freedom, here 
the gravitational ones. They are supposed to parameterize the reduced phase space and play 
a role of (partial) observables in quantum theory. It must be noted, however, that the com- 
mutation relations between these partial observables are postulated. They are parachuted 
from the kinematical phase space in an ad hoc manner. Since the observables are physical 
only in the constraint surface, any kinematical degree of freedom, Pj, forms the following 
equivalence class: 

Vi ~ V[ ^ V,^V[ ^ V, = Vl + Q (77) 

where Cj is a constraint. But then the Poisson bracket between the equivalence classes on 
the constraint surface is quite easily shown to be ill-defined (non- unique) : 

{V, + Q, V, + C,} = {Vi,V,} + {Q, V,} + m, C,} + {Q, C,} (78) 

where only the term {C,, Cj} ~ weakly vanishes. The terms {Cj, Vj} and {Vi, Cj} do not 
vanish in the constraint surface and make the Poisson structure ill-defined. 

The non-existence of the Poisson bracket in the constraint surface S was proved in the 
previous subsection and is due to the fact that the constraint surface is odd-dimensional. 
Therefore, in order to encode dynamics into quantum theory, one needs to postulate the 
Poisson bracket between partial observables or between a partial and Dirac one, so that 
relation fTMl) is fixed. Having this done, the time parameter is determined as the only partial 
observable which commutes with all the other, partial and Dirac, observables. Then, the time 
parameter t : S ^-^ R, together with the Dirac observables, introduces the 'no-Hamiltonian' 
parameterization of the constraint surface S = R x Vr. In the scalar field case t = <p, one 
has: 

{<l),V,}s = (79) 

Observe once again that in a constrained system, only the Poisson commutation between 
any pair of Dirac observables whatever its parametrization (see Appendix [B]) is well-defined: 

{V, + a, V, + C,} = V,] + {P„ C,} + {C„ V,] + {C„ C,] ^ {V,, V,} (80) 

since a Dirac observable by definition commutes weakly with a constraint, i.e. {Vi, Cj} ~ 
(in opposition to what happens in (!78|) ). 

To sum up, in the Dirac quantization, the problem of time persists. As before, this is 
so due to the fact that time- dependent observables are gauge-variant quantities and their 
Poisson commutation relations are undefined. In this procedure, an ambiguous Poisson 
structure can be postulated by parachuting the Poisson structure of 2r2 — 2 kinematical 
degrees of freedom of ambiguous choice. 



V. CHOICE OF TIME AND SPECTRA OF PARTIAL OBSERVABLES 

So far we have showed that the Poisson bracket between partial observables depends on 
the choice of time. In what follows, we will show how this fact affects spectral properties of 
time-dependent quantities. Let us study the Hubble observable in a fixed direction (1^ : 

H, = ^1^"^^ + '^^"'^ (81) 
V2(7 + COS0) tcos 
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First we notice that the evolution of the Hubble observable, in the cosmological time tcos, is 
not canonical. For from (IHTj) . we have: 

— + + (82) 



2 + 8 sin + 14 COS ' t, 
which after substituting for d(f) in pUj) leads to 



v/2//i(7 + cos0)^ 

= ■ ■ ■ + ^ , g ■ , , , . A dtcos (83) 

/ + 8 sm + 14 cos (p 

which according to the theory of canonical transformations (see the definition in the previous 
section) should have the form of ( 169|) . so that the one-form 

v/2//,(7 + co=*)= <iO. (84) 



2 + 8 sin + 14 cos ( 



should be equal to the derivative of a true Hamiltonian, cIHt, which here would be a gen- 
erator of the canonical motion of Hi. The existence of such a generator would allow us to 
construct a quantum theory in which the evolution of Hi would be unitary. This apparently 
does not hold, since the above one-form is not closed. In what follows, we will replace the 
cosmological time with a different clock variable. 

For brevity, let us restrict to G (0, and redefine time tcos in two ways: 

1 — cos c^) -I- 2 sin (/) / tt . , tt ^\ . . 

= /^rJ, IV^ IT sm t2 + O2 86 

V2(7 + cos0) \<Pi <Pi J 

Both ti and ^2 are well-defined, with ti E (— Oi, 00) for a given Oi and t2 E (— ^ ^) 

for a given Oi and (p. The times are half-lines and the left endpoints signal the singularity. 
Now, the Hubble observable flHTj) takes the form: 

^1 = ti + Oi = ^sin(^<^)t2 + 02 ^^^^ 

The difference between the above formulae is due to the different choice of time, i.e. a 
different parametrization of the constraint surface, though the function Hi : S ^ R itself 
remains unchanged, i.e. attaches real numbers to points in the constraint surface uniquely. 

Let us see that definitions f l85p and (jHSD lead to canonical motion in Hi. First let us 
examine the evolution in ti. Introduce: 

Ti:=ti + Oi (88) 

so that 

^1 = ^ (89) 



This phase space sector consists of all the solutions modulo the interchange of the axes in the homogenous 
leaf. 
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and 



Us = dTi Adr + d02 A dtp + dHr A dti 



(90) 



where Ht = r is a true Hamiltonian, which generates the canonical motion of Hi, in ti. For 
each constant time shce the symplectic form (1901) may be inverted to the Poisson structure: 



{-,■} 



d ■ d ■ d-d 



+ 



d ■ d ■ d ■ d ■ 



dTi dr dr dTi dO 



ti=const 

Now let us examine the evolution in t2- Introduce: 



To 



TT TT 

— sin(— 0)^2 + O2 
01 <Pi 



so that 



and 



Hi 



1 



Us = dOi Adr + dT2 A dtp + dHr A dt2 



(91) 



(92) 



(93) 



(94) 



where Ht = — cos(^0) is a true Hamiltonian, which generates the canonical motion of Hi, 
in t2- For each constant time slice the symplectic form flM|) may be inverted to the Poisson 
structure: 

d ■ d ■ d ■ d ■ ^ d ■ d ■ d ■ d ■ 

t2=const ^' ^' 9T2 dtp d(P dT2 



{-,■} 



dOi dr dr dOi dT2 



A. Warm-up: Quantization of Dirac observables 

We have identified the four-dimensional space of Dirac observables, (r, cp, Oi, O2) G -R+ x 
(0, (pi) X Rx R, equipped with the symplectic form ur identical with us in ( 140|) . By inverting 
the form in the reduced phase space one can find the Poisson bracket: 

{Oi,r} = l, {02,^} = 1 (96) 

We assign the following operators to the corresponding Dirac observables: 

r^r:=r, Oi^di:=idr, ^ := 0, O2 ^ 62 := id^ (97) 

so that the following algebra homomorphism is obtained: 

{0^} = -^[di,r], {0^} = -^[d2,^] (98) 
and the physical Hilbert space is defined as: 

-Hphys = L\R+ X (0, tPi),drdtP) (99) 

The formally self-adjoint operator id^, defined on a half-line, can be made into a self-adjoint 
operator by the following assignment: 

idr := (100) 
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where the Laplacian —d^ acts on the closure (in the 'operator norm') of the following space 
13: 

{^eUphys- ^(r = 0) =/i-9,^/;(r = 0)} (101) 

where ^ > enumerates unitarily inequivalent self-adjoint realizations of —d"^. The values 
= and fi = oo correspond to the Dirichlet and Neumann condition, respectively. The 
spectrum reads: 

SpR^{idr} = R+ (102) 

The formally self-adjoint operator id^i^ is also unbounded and enjoys many unitarily in- 
equivalent essentially self-adjoint realizations in the space 17l |: 

e Uphys : ^(0) = V^(0i)e*''^n (103) 

where and (pi are the boundary points of the domain. The parameter k G [0,1) enumerates 
unitarily inequvalent representations of id^i^. Its spectrum reads 

sp{id^} = ^ + ^n, n e Z (104) 

<Pi (Pi 

This completes the quantization of all the gauge-invariant observables in the Kasner 
universe. Since the notion of evolution is absent in this model, one cannot ask questions 
about the fate of singularity. 



B. Quantization of geometrical observables in two different time variables 

1. Case ti 

Let us use the Schrodinger representation to study the Hubble observable defined in ( 189|) 
with the Poisson structure given in (!9T|) : 

r ^ f := -i-^, Ti ^ Ti := r, := 0, O2 ^ 62 := -i-^ (105) 

dr dcp 

which act on the Hilbert space 

Hphys = L\R+ X (0, 0i), drd(t)) (106) 
so that the following is satisfied: 

{TC^} = i[ri,f], {a^} = ^-[d2A] (107) 

The true Hamiltonian reads: 

= -zA (108) 
dr 

The operator —i-^ was discussed in the previous subsection, and we showed that it is essen- 
tially self-adjoint in the domain fllOip . Thus, due to the Stone- von Neumann theorem there 
exists a unitary operator 

U = e-''^^ (109) 
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which proves that the dynamics of the system is well defined for all t G R. This means that 
the singularity which is reached in finite time in classical theory is resolved at quantum level. 
Note that the quantum operators associated with physical measurements (like connection 
or curvature) may be unbounded. These circumstances, however, are common in quantum 
theory and do not spoil the singularity resolution. 
Let us move to the quantization of Hi: 

Hi = ^ = -, reR+ (110) 
Ti r 

which is a well-defined self-adjoint (unbounded) operator on a half-line with the continuous 
spectrum sptA^i) — 



2. Case t2 

Let us use the Schrodinger representation to study the Hubble observable defined in ( p3|) 
with the Poisson structure given in (195|) : 

r ^ f := -i^, Oi^di:=r, 0^0 := (j), ^ := -z-^ (111) 
dr dcp 

which act on the Hilbert space 

U^hys = L\R+ X (0, 0i), drd(P) (112) 
so that the following is satisfied: 

{0^} = -[di,r], {T^} = -[f,,^] (113) 

The true Hamiltonian reads: 

Ht = - cos (—(/)) (114) 
01 

The operator — cos(^0) is bounded, symmetric and hence self-adjoint in the Hilbert space. 
Due to the Stone-von Neumann theorem there exists a unitary operator 

U = e-'^^^ (115) 

which proves that the dynamics of the system is well defined for all t & R. This again means 
that the singularity which is reached in finite time in classical theory is resolved at quantum 
level. 

Let us move to the quantization of Hi: 

Hi = ^ (116) 

The operator: 

^ d , ^ 

= -'T^ ("^) 

has been already discussed and is a well-defined self-adjoint (unbounded) operator with the 
discrete spectrum sp{id^} = ^ + ^n. Using the spectral theorem and requiring positivity 

of Hi, we find that the spectrum of the Hubble observable is discrete, bounded and reads 



spt2{Hi 



neZ 
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VI. CONCLUSIONS 



Motivated by the BKL scenario, we started this paper aiming at deriving a quantum 
theory of the Kasner epoch. This task seemed to be feasible with the use of the reduced 
phase space method: we defined the kinematical phase space, in which we solved Hamilton's 
equation of motion; then we identified Dirac observables and their algebra; finally, we arrived 
at the physical solutions of the classical theory in terms of Dirac observables and a clock 
variable. At this point, the quantization is usually performed. However, we realized that the 
choice of clock variable determined the functional dependance of time-dependent quantities 
on Dirac observables. In the rest of the paper, we studied the consequences of this fact. 

We showed that in addition to the usual ambiguities of quantum theory, in a Hamiltonian 
constraint system like general relativity, there is also another ambiguity related to the choice 
of the clock. We have managed to clarify the procedure of encoding evolution in quantum 
gravity. It turned out that the procedure could be identified with inverting a singular matrix. 
The matrix is the induced two-form on the constraint surface a;^ and its inverse, in a sense, 
still exists but is ambiguous and represents the Poisson bracket. Different choices of the 
inverse lead to canonically inequivalent classical theories, which are then quantized. Then 
the dependance of quantum physics on the choice of time variable can be stated as follows: 
The Poisson bracket associates a canonical transformation with each phase space function. 
In classical theory, the Poisson bracket is an auxiliary mathematical structure, which does 
not affect the physical (observable) content of the theory in the sense that it does not matter 
which canonical transformation corresponds to which phase space function as long as the 
equations of motion are equivalent. In quantum theory, however, the physical quantities are 
associated with operators on the Hilbert space, in which they act as generators of unitary 
transformations. The correspondence between classical and quantum theory is the one be- 
tween the canonical transformations and the unitary ones. Therefore, the Poisson structure 
is not auxiliary but essential for this correspondence and thus all ambiguities in canonical 
formulations of classical theory are expected to lead to quantum theories with different phys- 
ical content. Thus, for each choice of time, one may find a distinct spectrum for a given 
partial observable. We have shown that this is the case for a directional Hubble parameter 
for which we were able to obtain both continuous and discrete spectrum, depending on the 
choice of the clock variable. Which spectrum, if any, is the correct one? Answering this 
question is beyond the framework of general relativity and usual quantum mechanics. 

In the view of this result, one should construct all the quantum theories treating all the 
possible time parameters on equal footing. The relation between the resultant canonically 
inequivalent theories surely deserves further study. This will be the subject of next papers 
by the present author. 

What is the physical interpretation of the obtained result that theories of gravity with 
different time parameters are canonically inequivalent? Is the evolution of the universe not 
an objective reality but merely an impression perceived by an observer? Is it possible for 
every observer to 'measure' his own story of the big bang? At this moment, it is very 
tempting to speculate that perhaps the fundamental quantum theory we all are looking for 
is not a gauge theory - the gauge invariance is obtained only in the classical limit, in which 
the uncertainty principle is removed. This point of view is strongly supported by the fact 
that one can quite easily reformulate a classical non-relativistic system in a gauge-invariant 
manner (see e.g. [lo[). The natural question then arises: among all the partial observables, 
which one is always classical? 
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Appendix A: Calculation of connection and curvature 

We consider a manifold Ai equipped with the metric: 

i 

where i = 1,2, 3, and on which we introduce the vector fields and the dual 1-forms: 

Ndt , cr* = ttidx^ 
so that the following relations hold: 

where /i, = 0, 1, 2, 3. We assume that fisl ]: 



eo = -T^Ot , Ci = —di, a 
N a. 



(Al) 



(A2) 



(A3) 



(A4) 



and obtain the connection Q: 

da = aQ , Q 

which is a matrix of 1-forms. It is uniquely determined from the condition u^^ + w 
and equals: 

/ . \ 



-^a^ 

_a3_ 3 















/ 



From Vt we compute the curvature matrix B: 

e = w,\\ 



dVL-^f 



which is equal to 

The Riemann curvature tensor is now given by: 



e 



Na:i \ N ) t \ 

^ Nai > \ Nas " 
/ji^Wji^N 2 3 





(A5) 
= 0, 

(A6) 



(A7) 



(A8) 



(A9) 
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from which we calculate the Ricci tensor R^/s: 



/ Roo \ 



7?22 

\ R33 J 



(AlO) 



where i?oo = - ]^ (|) , and i?,, = ^ (| ) (A) i^) ' ^hus, the Ricci scalar: 



(All) 



In the Friedman cosmology, it is common to use the Hubble parameter H and the deceler- 
ation parameter q, which in the case of Kasner universe can also be introduced, for each of 
the three directions separately: 



1 (a^\ 
Na, \ N J , 



m 



(A12) 



Note that these parameters determine all the components of the curvature matrix 6. 



Appendix B: Parametrization of Dirac observables 

It is important to realize the precise notion of Dirac observable. It is sometimes claimed 
that Dirac observable is a kinematical phase space function, which commutes weakly with 
the Hamiltonian constraint. This is not precise enough. In a constrained system, the physical 
motion is realized in the constraint surface, let us say H = 0. This equality implicates the 
embedding of the constraint surface into the kinematical phase space, E : S ^ V. One 
cannot pull back functions from S to functions on T', since the range of E is restricted only 
to H = This means that the Dirac observable cannot be defined as a function of the 
kinematical phase coordinates. The proper definition is as follows: Dirac's observable Vi 
is a constraint surface function, which commutes with the Hamiltonian constraint H, i.e. 
XH{Vj)\^ = 0, where Xh satisfies uj{-^Xh) = dH. 

Let us find functions Od, which satisfy: 



{Od,H} = 



(Bl) 



In the Kasner universe the above equality in the convenient gauge N — A\/2\/ X1X2X3 takes 
the following form: 



dOp 
dXi 







and we obtain: 
where 



Qij = (r^ + r,)(ri + r^) in IP^I - (r^ + r,)(r,- + r^) in \Pj\ 



(B2) 

(B3) 
(B4) 
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and satisfy the relation: 

^ij — ^kj = ^ik (B5) 

These are solutions to fIBip in a given gauge. Their restriction to the constraint surface 
f|T8|) is gauge-invariant and gives the four-dimensional space of Dirac observables, Vr. We 
could alternatively define the Dirac observable as equivalence classes with the equivalence 
relation: 

~ ^ Od^O'^ ^ Od = 0'^ + C (B6) 

where '~' denotes 'equals on the constraint surface' and C ~ is a constraint. This, how- 
ever, leads to the ambiguity in 'parameterization' of Dirac observable as there is no natural 
projection from the kinematical phase space to the reduced phase space. The important 
issue here is that the Poisson bracket is well defined between the equivalence classes (see 
equation flHOl) ). 

Appendix C: The Kasner metric in terms of the cosmological time and Dirac ob- 
servables 

Let us express the physical solutions in terms of Dirac observables and cosmological time 

tcos- 

2 2 — 2 cos 0+4 sin 

-2^ r(7 + cos0) \ 7+^Hi0 



ds = -dt^^^ + 2a^[ ^^-j= — teas ) (dx 



r(7 + cos 0) 



2 — 2 cos (p — A: sin q 



2 3+5 cos (p 

_9/ r(7 + cos0) \ ■ 7+cos 



where 



{Oj (7 sin 0-4 cos 0) + i 02(1 + 3 cos 0+4 sin 0) } {-5-9 cos 0-4 sin 0} 
g 2(7+cos 0) (5+3 COS 0— 4 sin 0) 

Ox (^ sin + 10 cos 0+^ sin cos 0— 14sin2 + 6 cos^ 0)+ i O2 ( 15- 12 sin 0+cos 0-12 sin cos 0- § sin2 0) 
0,2 = e (7+cos 0)(5 + 3cos0-4sin0) (C2) 

Oi(-20 cos 0+16 sin cos 0-12 cos2 0)+ 1 O2 (- 16+20 sin 0+8 cos + 12 sin cos 0+16 cos2 0) 
^ g (7+COS 0)(5+3 COS — 4 sin 0) 

The minus sign in front of the time, —tcos, comes from our convention that as time grows 
the universe approaches the singularity, that is tcos G (— oo,0). Now, we can relate the 
parameters pi, and p^ occurring in the metric ([5]) to the Dirac observables: 

Fi 2 — 2 cos (p + 4 sin (p 
7 + cos (p 

T2 2 - 2cos(;/) - 4sin0 

P2 = ^^-fr = ^ 7 (C3 

Z^i 7 + cos 

Fa 3 + 5 cos (j) 
Y^i^i 7 + cos 

We note that all the Dirac observables play a role provided that the topology of the 
universe is compact. If, for instance, we set S infinite in all directions, then cj) is the only 
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Dirac observable, since the values of the scale factors are non-physical, and consequently 
r, Oi and O2 are not Dirac observables any longer. 
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